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Abstract 

We construct a low-energy effective Lagrangian describing zero-temperature supersolids. 
Galilean invariance imposes strict constraints on the form of the effective Lagrangian. We identify a 
topological term in the Lagrangian that couples superfluid and crystalline modes. For small super- 
fluid fractions this interaction term is dominant in problems involving defects. As an illustration, 
we compute the differential cross section of scatterings of low-energy transverse elastic phonons by 
a superfluid vortex. The result is model-independent. 
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Introduction. — The possibility of superfluid behavior in sohds was considered by Andreev 
and Lifshitz (AL) in a seminal paper of 1969 [1]. They proposed that point defects in 
a ^He crystal may become quantum at low temperatures and form a Bose condensate. 
The possibility of superfluid solids ( "supersolids" ) was also conjectured in Refs. [2, 3]. For 
a long time, the experimental search for a supersolid phase of helium was unsuccessful. 
Recently, however, Kim and Chan claimed that superfluid behavior may have been observed 
in solid ^He [4, 5]. This has stimulated renewed interest in the supersolid phase. The 
superfluid fraction was found in Ref. [5] to be of order 10~^, considerably larger than what 
was theoretically expected (< 10~^ [3]). 

In the following, we assume that the supersolid phase does indeed occur. We shall not 
discuss the microscopic origin of such a behavior. Rather, we concentrate our attention 
on the theoretical description of the low-energy dynamics of the supersolid phase. One can 
expect that, regardless of the details of the microscopic mechanism underlying the supersolid 
state, its low-energy dynamics is simple and universal, and is dictated by symmetry principles 
alone. This is because the low-energy degrees of freedom of supersohds are the Nambu- 
Goldstone bosons arising from spontaneous breaking of translational symmetry and the 
U(l) symmetry generated by the conserved particle number. In fact, Andreev and Lifshitz 
have constructed the hydrodynamic equations describing a supersolid (in the limit of small 
strain), based on conservation laws and Galilean invariance only. 

In this Letter, we restrict ourselves to zero temperature and derive an alternative de- 
scription of supersolids based on a low-energy effective Lagrangian [6]. Such a description 
is possible since dissipative effects disappear at T = 0. The effective Lagrangian description 
holds several advantages over the one based on the hydrodynamic equations. As it is gener- 
ally true, the Lagrangian formulation enables straightforward application of field-theoretical 
techniques such as Feynman diagrams. In our particular case, the Lagrangian also elucidates 
the appearance of a certain "topological" interaction term, which is important in processes 
involving defects. 

The form of the effective Lagrangian is constrained by various symmetries, among which 
an important role is played by the Galilean invariance. For liquid supcrfluids, the most 
general Galilean-invariant effective Lagrangian was constructed by Grcitcr. Wilczck, and 
Witten [7]. For supersolids the possible structures appearing in the Lagrangian is richer 
than in supersolid ^He, but are still rather restrictive. 
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Degrees of freedom. — The fields appearing in the effective Lagrangian are the four Gold- 
stone bosons which appear due to the spontaneous breaking of the U(l) particle number 
symmetry and the translational symmetry along three spatial directions. One of these fields 
is therefore the phase of the superfluid condensate, 6; under the action of the particle num- 
ber, it transforms as 6* — > 6* + a. 

The three remaining fields are translation- breaking scalars X", a — 1,2,3, which can be 
introduced as follows [8]. Imagine a system of coordinates X", a — 1, 2, 3 which is frozen in 
the body of the sohd. (The choice of the system is completely arbitrary, but for simplicity 
one can choose it to coincide with the coordinates in our reference frame x', i — 1, 2, 3 if 
the solid is at equilibrium at some arbitrarily chosen external pressure Pq.) When the solid 
moves, this system of coordinates also moves, so if one follows one particular material point 
in the solid, its coordinates in the X system remain constant. In general, the coordinate 
system X is curved. The time history of the solid is completely characterized by three 
functions X°'{t, x), which give the coordinates, in the comoving frame, of the material point 
that is located at the position x at time t. X"(t,x) are the fields that enter the effective 
theory together with the U(l) phase 9. 

One can expand the fields around the ground state 

e^Hot- ip, X^^x"- -u", (1) 

where (p and fluctuate around zero. Here /ig is the chemical potential at pressure Pq] 
is the usual displacement vector. The superfluid velocity and the strain of the crystal are 
related to the flrst spatial derivatives of 9 and X"-: 

V, = - V(^, diX'^ = - diu'' (2) 
m 

(in most of this Letter, Yi—V). 

The density of lattice sites is a constant in X space. We denote this constant by no. In 
the X space, the density of lattice sites is 

no det %X-\ = '^e,i^e^,d,X'^diX'>d^X\ (3) 

For an ordinary crystal, nonsuperfluid and without defects, this coincides with the particle 
number density. 

Derivative expansion. — The effective Lagrangian should be invariant under the U(l) par- 
ticle number symmetry, 9 ^ 9 + a, and coordinate shift in the frozen frame, X" X'^ + a"'. 
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Thus the Lagrangian contains only time and coordinate derivatives of 9 and X", but not 9 
or X"- by themselves. 

In order to discuss the low-energy regime, we follow the standard effective field theory 
philosophy and perform a derivative expansion of the effective Lagrangian. Because 9 and 
X"- do not appear without derivatives, there are two different ways to perform the derivative 
expansion. One possibility is to assume that 9 and X"- are slowly varying functions of space 
and time. Alternatively, one can assume that the first (temporal and spatial) derivatives of 
9 and X"- vary slowly. Clearly, the second alternative is more general, since it allows for the 
first derivatives of 9 and X"- to be large. In particular, the superfiuid velocity and the strain 
of the crystal do not have to be small. In our subsequent discussion, we will therefore assume 
that 9, di9, X" and diX"- may be not small, but vary slowly in space and time over distances 
large compared to all microscopic length scales, such as the superfiuid healing length. 

Keeping only leading-order terms in the derivative expansion, the Lagrangian is a function 
of the first derivatives of fields, 

jC^jC{9,di9,X'',diX^). (4) 

This Lagrangian, in general, contains terms to all orders of fields. In each term in the series 
expansion over fields 9 and X"', one keeps the lowest possible number of derivatives equal to 
the number of fields. 

Rotational invariance. — The effective Lagrangian should be invariant under spatial ro- 
tations. The fields X", despite being a three-component field, transform under spatial 
rotations hke scalar fields. This is because X" are the coordinates of the internal system 
frozen in the solid body, which are not rotated with the axes of spatial coordinates x*. There- 
fore the Lagrangian (4) is a function of the following rotationally invariant combinations of 
arguments, 

£ = £{9, X\ {di9f, di9d,X\ u''^) , (5) 
where we introduce the notation 

u"^ = diX'^diX^ = 5"'' - ay - d^u" + diu'^diuK (6) 

It can be shown that any rotationally invariant function of first derivatives of 9 and X" can 
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be written as a function of the parameters staying in Eq. (5). In particular, 

eijkdiX^'djX^dkX^ = e^'^VdetM , (7) 
eijkdi0djX''dkX^ = e''^^Vd^uJdi0diX'^, (8) 

where det u is the determinant of the 3x3 matrix u"'', and u~} is the inverse matrix of u: 



ah 

Kb'^bc = Sac- 
Galilean invariance. — Further constraints on the effective Lagrangian follow from 
Galilean invariance. In a nonrelativistic theory where all particles have the same mass, 
the momentum density is proportional to the particle number flux, 

T'' = mj\ (9) 

where m is the mass of the '^He atom. The momentum density and the particle number flux 
are found from the Lagrangian by using Noether's theorem, 

T^i ^ -^diO - ^d,X\ f = . (10) 

09 dX^ d{di9) ^ ' 

The most general form of the Lagrangian which is consistent with rotational [Eq. (5)] and 
Gahlean invariance [Eq. (9)] is 

>c = i:(/i,w",M"''). (11) 

Here u"'^ was defined in Eq. (6), and ^ and w"" are 

H = e- ^^f^ , w" = -X" + -diOdiX^ (12) 
2m m 

The variable /x appears in the Lagrangian treatment of superfluids [7]: the crystal structure 
is absent there and the Lagrangian is a function of alone. Physically, /x is the local chemical 
potential as measured in the frame moving with the superfluid velocity. The meaning of 
can be made clear by expanding it, 

w = ix+(v, • V)u-v,. (13) 

At the linearized level w"' is the difference between the velocity of motion of the crystal 
lattice and the superfluid velocity. 

Connection to the AL hydrodynamic theory. — The zero-temperature AL hydrodynamic 
equations can be derived from the Lagrangian (11) as the equations of conservation of particle 



number, energy, and momentum, written using a particular set of variables. We introduce 
the total density p, the superfluid velocity Vs, and the vector which corresponds to the 
normal velocity in the AL theory, 

P = rn— , Vsi = —Oiif, Vni = -w^u . (14) 

Actually v„ is the velocity of the crystal lattice; at the linearized level v„ = u. At zero tem- 
perature we do not have a normal component distinct from the crystal lattice. Furthermore, 
we introduce the momentum density p and the energy density e in the frame Vg = 0, 

(1^) 

and an auxiliary tensor Aja, 

dC dC 
-^- = 2^^.^' + K.-^..)^. (16) 

After some algebra, we find that the momentum density can be written as 

T^^rnf^pvsi+Pi, (17) 

the energy density and the energy flux as ^ 

= ^ + p-v. + e, (18) 

and the differential of e as 

de = ^dp + (v„ - V,) • dp + Kadidiu'^) . (20) 
m 

Equations (17)-(20) are identical to the corresponding AL equations at T = 0. The stress 
tensor can be transformed into the form 

= pVsiVsk + VskPi + VniPk + 5ik II— + (V„ - V^) • p - £ - A^fc + XiadkU"" , (21) 

L m J 

which almost coincides with the corresponding AL expression. The only difference is the 
last term on the right hand side of Eq. (21), which is nonlinear in strain and was neglected 



^ Note that T"* ^ T**^, since in our nonrelativistic theory energy does not include rest mass. Our normal- 
ization of the chemical potential differs from that of AL by a factor of m. 



6 



in Ref. [1]. Moreover, from the definition of Xia, Eq. (16). one can derive the following 
relation: 

Aifc - Xki - XiadkU"- + XkadiU"- = {Vni " Vsi)pk " {Vnk " Vsk)Pi , (22) 

which coincides, up to the two terms nonlinear in strain on the left hand side, to an equation 
postulated in Ref. [1] for T*'^ to be a symmetric tensor. 

As one can see, the effective Lagrangian provides an extremely compact encoding of the 
hydrodynamic equations. Terms nonlinear in strain which were omitted in the AL theory 
are fully kept in the Lagrangian. Moreover, for quantum problems (such as scatterings of 
phonons) it is easier to work with the Lagrangian than with the field equations. We now 
show that the Lagrangian contains a special topological term that is important for scattering 
off defects. For this end, we first discuss the nonsuperfluid hmit of the Lagrangian (11). 

The nonsuperfluid limit — The superfiuid fraction ps/p of solid helium-4, if nonzero, is 
much smaller than one: experiments [4, 5] indicate a value of order 10~^, while theoretical 
arguments [3] suggest Ps/p < 10~^. We shall, therefore, concentrate on the hmit Ps <^ p. 
To start, let us discuss the limit of vanishing superfiuid density. 

One expects that the nonsuperfluid crystalline state is realized as some particular limit 
of the supersolid state, where 9 decouples from the dynamics of fields. This might seem 
nontrivial, since in Eq. (11) the time derivative of X" enters the Lagrangian in the combi- 
nation which involves 9 [Eq. (12)]. However, it is possible to achieve such a decoupling. 
Consider the following Lagrangian: 

c = poVd^ ilKl^'y^' + -)- ^(«"') , (23) 

where po is some constant with the dimension of mass density, p and w are defined in 
Eq. (12), and V is an arbitrary function of the strain u""^ consistent with lattice symmetry. 
Using the identities (7) and (8), this Lagrangian can be transformed into the form 

£ = ^Vd^u-lJC^X' - V{u''') + ^e>'''^''eabcd^9d,X''dxX%X''. (24) 

The Greek indices in the last term are spacetime indices which run over t, x, y, z; e'^"^^ is the 
completely antisymmetric tensor defined so that e°^^^ — 1. The phase 9 appears only in the 
last term of the Lagrangian, which is a full derivative; thus 9 decouples from the dynamics. 
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The last term, which will be called the "topological term," still plays a useful role: the 
particle number current computed from (24) by using Noether's theorem arises completely 
from this term, 

f = ^e^''^''eabcd.X'^d,X%X^, (25) 

from which we see that po is the total mass density in equilibrium where X"' — x"'. The 
current (25) is trivially conserved. Note that the first term in Eq. (24) can be written as 
mfj^l {2j'^), which is what one expects for the kinetic energy from Galilean invariance. This 
fact shows that the construction (23) is unique. 

Now if one allows the superfiuid fraction ps to be nonzero and small, then there are 
additional terms proportional to Ps added to the Lagrangian, which makes 9 a dynamical 
field. The topological term continues to be present in the Lagrangian with a coefficient 
which differs only slightly from p/(6m). This term is responsible for low-energy scattering 
of elastic waves by a superfiuid vortex, as we shall see. 

Scattering of elastic waves by a superfiuid vortex. — Let us now use the effective Lagrangian 
to compute the scattering cross section of elastic waves by a superfiuid vortex. In the 
presence of a vortex (/? is a multivalued function, and ~ V</7 is singular at the vortex core. 
Due to the multivalued nature of 9, the topological term is no longer integrated by part to 
zero. Integrating by part, this term can be written as 

Aop = -£^e''''^''e'''%d,ipX''d^X%X'. (26) 

Expanded over small perturbations, this expression contains a term proportional to u • v^, 
which was identified in Ref. [9] in the context of an Abrikosov vortex in a crystal. For 
definiteness, consider a vortex located at x — y — and stretched along the z direction. 
The field of the vortex has 

dxdy(p - dyd:c^ = 2TiS{x)S{y) . (27) 

The topological term is localized on the vortex core and has the form 

jC = -7r—5{x)5{y){uxUy - Uyiix) + 0{u^) . (28) 

Note that the leading term in this Lagrangian contains two powers of u but only one deriva- 
tive, so it cannot be canceled by an unknown interaction of the vortex core with the elastic 
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FIG. 1: The scattering of transverse elastic phonons by a superfluid vortex. 



waves. Moreover, other terms that couple 6 and X" are expected to be proportional to ps 
and are negligible. 

The process that will be considered has the following kinematics (Fig. 1). A transverse 
phonon with momentum fik and linear polarization e falls perpendicularly onto a vortex. 
We are interested in the probability of its scattering into the state with momentum fik' and 
polarization e' . For simplicity, we assume that both e and e' lie in the plane perpendicular 

to the vortex and that the solid is isotropic. The matrix element of the process is 

M = — (exe')-z, (29) 

m 

from which we find the differential cross section per unit vortex length, 

ksin^e, (30) 



dedl 2m?vl 

where v±_ is the speed of transverse elastic waves. This result is a model-independent pre- 
diction of the effective Lagrangian approach, valid at small k and small superfluid fraction. 

Qualitatively, the differential cross section has a linear dependence on k and is maximum 
when the scattering angle 9 is 90°. 

Conclusion. — We have found the most general effective Lagrangian describing the low- 
energy dynamics of supersolids. We show that, in the limit of small superfluid density, the 
Lagrangian contains a topological term which has a fixed coefficient. From this term we 
computed the cross section of scattering of transverse phonons off a superfluid vortex. If the 
supersolid state is realized in '^He, this prediction is, in principle, verifiable. 

The formalism used in this paper can be extended to relativistic systems, e.g., for describ- 
ing the crystalline superfluid phases of quark matter [10]. Instead of Galilean invariance, 
one requires relativistic invariance of a theory of four Goldstone bosons 9 and X". Such a 



theory, when coupled to gravity, gives rise to a gravitational analog of the Anderson-Higgs 
mechanism. Indeed, theories of this type have been proposed recently as an infrared modi- 
fication of gravity that gives the graviton Lorentz- breaking mass terms [11, 12]. Prom this 
point of view, the modification of gravity considered in Refs. [11, 12] can be interpreted as 
an effect coming from a supcrsolid dark matter sector. 

The author thanks A. Andreev, J. Erlich, P. Kovtun, M. A. Stephanov, and D.J. Thouless 
for discussions. This work was supported by DOE grant DE-PG02-00ER41132 and the Alfred 
P. Sloan Poundation. 
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